Goldstone bosons in presence of charge density 

Tomas Braune 
Department of Theoretical Physics, Nuclear Physics 
Institute ASCR, 25068 Rez, Czech Republic 

Abstract 

We investigate spontaneous symmetry breaking in Lorentz-noninvariant theories. Our general 
discussion includes relativistic systems at finite density as well as intrinsically nonrelativistic sys- 
tems. The main result of the paper is a direct proof that nonzero density of a non-Abelian charge 
in the ground state implies the existence of a Goldstone boson with nonlinear (typically, quadratic) 
dispersion law. We show that the Goldstone boson dispersion relation may in general be extracted 
from the current transition amplitude and demonstrate on examples from recent literature, how 
the calculation of the dispersion relation is utilized by this method. After then, we use the general 
results to analyze the nonrelativistic degenerate Fermi gas of four fermion species. Due to its 
internal SU(4) symmetry, this system provides an analog to relativistic two-color quantum chro- 
modynamics with two quark flavors. In the end, we extend our results to pseudo-Goldstone bosons 
of an explicitly broken symmetry. 
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I. INTRODUCTION 



The basic result concerning the physics of spontaneous breaking of global internal con- 
tinuous symmetries in Lorentz-noninvariant systems was achieved by Nielsen and Chadha 
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thirty years ago They showed that, under certain technical assumptions (the most im- 
portant one being that of a sufficiently fast decrease of correlations, which ensures that no 
long-range forces are present), the energy of the Goldstone boson (GB) of the spontaneously 
broken symmetry is proportional to some power of momentum in the low-momentum limit, 
and classified the GBs as type-I and type-II according to whether this power is odd or even, 
respectively. Their general counting rule states that the number of type-I GBs plus twice the 
number of type-II GBs is greater than or equal to the number of broken generators. 

Two classic examples of systems where the type-II GBs occur are the ferromagnet and 
the so-called A phase of superfluid 3 He. The issue of GBs in Lorentz-noninvariant theories 
has regained considerable interest in the past decade due to the discovery of several other 
systems exhibiting the existence of type-II GBs. These fall into two wide classes: Relativistic 
systems at finite density and intrinsically nonrelativistic systems. The first class includes 
various color-superconducting phases of dense quark matter — the color-flavor-locked (CFL) 
phase with a kaon condensate [2J, y] and the two-flavor color superconductor |4] — as well as 
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the Bose-Einstein condensation of relativistic Bose gases [5|, [6( and the relativistic nuclear 
ferromagnet |7j. The second class covers the before mentioned ferromagnet and superfluid 
helium, and the rapidly developing field of the Bose-Einstein condensation of dilute atomic 
gases with several internal degrees of freedom 8]. 

Besides the investigation of particular systems, a few new general results also appeared. 
Leutwyler IJ] showed within the low-energy effective field theory framework that the non- 
linearity of the GB dispersion is tightly connected to the fact that some of the conserved 
charges develop nonzero density in the ground state. As a matter of fact, the nonzero charge 
density induces breaking of time-reversal invariance, which in turn leads to the appearance 
of a term in the effective Lagrangian with a single time derivative. This term is then respon- 
sible for the modification of the GB dispersion. In addition to Leutwyler's results, Schaefer 
et al. 3| proved that nonzero ground-state density of a commutator of two broken charges 
is a necessary condition for an abnormal number of GBs. 

In our previous work 10] we showed that, at least within a particular class of systems 
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described by the relativistic linear sigma model at finite density, and at tree level, the 
converse to the theorem of Schaefer et al. is also true: Nonzero density of a commutator 
of two broken charges implies the existence of a single GB which has quadratic dispersion 
relation, i.e., is type-II. Moreover, there is a one-to-one correspondence between the pairs 
of broken generators whose commutator has nonzero density, and the type-II GBs. In Ref. 

ll| we analyzed the particular model of Miransky and Shovkovy [2|] and Schaefer et al. 

3J at one-loop level and asserted that the one-loop corrections do not alter the qualitative 
conclusions achieved at the tree level. 

The purpose of this paper is to extend the previous results in a model-independent way. 
We are going to present a general argument, based on a Ward identity for the broken 
symmetry, that nonzero ground-state density of a non-Abelian charge leads to the existence 
of a GB with nonlinear (typically, quadratic) dispersion law. The plan of the paper is 
following. In the next section we give a detailed derivation of the main result. In the 
subsequent sections we then apply it to several systems discussed in recent literature, and 
show how it facilitates the calculation of the GB dispersion relations. Next we analyze in 
detail Cooper pairing in a nonrelativistic gas of four degenerate fermion species. Finally, we 
show how the proposed method may be generalized to include explicit symmetry breaking. 

II. GENERAL ARGUMENT 



In Ref. [10j we provided a general argument that nonzero density of a non-Abelian 
charge leads to modified GB counting. In particular, we invoked the standard proof of the 
Goldstone theorem [TJ] to show that 

iUMfc^M = (0|[j»,g 6 ]|0) = 2.Im^(0|j a (0)|n)(n| J& °(0)|0), (1) 



see Eq. (1) in Ref. llOj and the discussion below it. Here f a bc is the set of structure 
constants of the symmetry group and the index n runs over different zero-momentum states 
in the spectrum. The vacuum expectation value (0|j°|0) represents the order parameter for 
symmetry breaking, while both currents j° and j® serve as interpolating fields for the GBs. 
Qb is the charge operator associated with the current j£. Note that whereas in Lorentz- 
invariant systems there is a one-to-one correspondence between the GBs and the broken 
currents, at nonzero charge density the situation changes. As Eq. (CQ) clearly exhibits, a 
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single GB now couples to the two broken currents whose commutator has nonzero vacuum 
expectation value. Before we expand this simple observation into a more detailed argument, 
we discuss what the involved current transition amplitudes tell us about the GB dispersion 
relation. 



A. GB dispersion from transition amplitude 

The matrix element which couples the GB state |7r(k)) to the broken current operator 
ia (0), may be generally parameterized as 

(0|^|7r(k))=z^F a (|k|)+^°G o (|k|). (2) 

By the second term, G a , we explicitly account for the possibility of Lorentz violation. How- 
ever, we still assume that the rotational invariance remains intact. The current conservation 



then implies the following equation for the GB dispersion (see also Ref. [13]), u;(k), 
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k 2 )F a + ^G a = 0. (3) 



Let us now explicitly mention a few special cases, distinguished by the limiting behavior 
of the functions F a and G a at low momentum. First, when G a = for all values of k, we 
recover the standard Lorentz-covariant parameterization of the transition amplitude (j2J), and 
the standard Lorentz-invariant dispersion relation, u = |k|. Second, suppose that G a /F a is 
small in the limit |k| — > 0, of order 0(|k|). Then the dispersion relation comes out linear, 
only the phase velocity is no more equal to one due to the Lorentz- violating term G a . Third 
and most importantly for the rest of the paper, when G a has a finite limit, lim^^o G a ^ 0, 
the dispersion relation takes the low-momentum form 

w(k) = k 2 ^. (4) 

Eq. (jl]) implies that the GB dispersion relation is nonlinear. However, in order to determine 
the corresponding power of momentum and thus classify the GB according to the Nielsen- 
Chadha scheme, we would have to further know the limiting behavior of F a . 

One should keep in mind that so far, we have provided just a parameterization of the 
GB dispersion relation in terms of the current transition amplitude, which only assumed 
rotational invariance and current conservation. In order to decide which of the possibilities 
outlined above is actually realized, we need more information on the dynamics of the system, 
i.e., on the functions F a , G a . 



B. Transition amplitude from Ward identity 

We now get back to Eq. (pQ). In order to turn it into a quantitative prediction of the GB 
dispersion relation, we replace the commutator with the time ordering, that is, consider the 
correlator of two currents, (Q\T{j%(x)j%(y)}\0) . Taking a divergence and using the current 
conservation, we arrive at a simple Ward identity, 



Of course, with the assumed rotational invariance the right hand side of Eq. may only be 
nonzero for v = so that it exactly encompasses our order parameter — the charge density. 

Throughout the paper we assume that the symmetry is not anomalous so that, in par- 
ticular, the naive Ward identity (jSJ) holds without any corrections due to, e.g., Schwinger 
terms. This assumption is partially justified by the models studied in the following sections, 
where our results prove equivalent to those obtained with different methods. In general, 
however, the validity of Eq. (jSJ) should be checked case by case. 

Next we use the Kallen-Lehmann representation for the current-current correlator in the 
momentum space, plugging in the parameterization (j2J). We find that the condition (j3J), as 
expected, ensures cancelation of the one-particle poles in the correlator. Taking the limit 
|k| — > 0, it follows that in order to get a finite nonzero limit as the right hand side of the 
Ward identity (jSJ) requires, the functions G a , Gb have to have a nonzero limit as well. The 
residua of the poles then yield the density rule G a Gl — = if a b c (0\jc\fy, or > 



where the Gs now stand for the zero-momentum limits of the respective functions defined 
by Eq. (j2J). Eq. (jSJ) together with Eq. (jlj) constitute the main result of the paper. They 
say that nonzero density of a non-Abelian charge in the ground state requires the existence 
of one GB with nonlinear dispersion relation. In general (in fact, in all cases the author is 
aware of), this dispersion is quadratic unless F a vanishes at k = 0. We will speculate on 
this possibility in the conclusions. 

A few other comments are in order here. First, the fact that the function G a has nonzero 
zero-momentum limit essentially means that the broken charge operator creates the GB even 
at zero momentum with a finite probability. This is in contrast to the Lorentz-invariant 




(5) 



lmG a G* b = ±f abc (0\j M, 



(6) 



5 



case. Second, in dense relativistic systems the function G a must be proportional to the only 
Lorentz- violating parameter in the theory, the chemical potential. An example will be given 
in the next section. Last, in view of Eqs. (jl]) and (161). it is tempting to conclude that the 
energy of the GB is inversely proportional to the charge density. This indeed happens in the 
ferromagnet [9(, but need not be true generally. Again in the next section, we shall see that 
both F a and G a may be proportional to the charge density so that it drops in the ratio. 



III. LINEAR SIGMA MODEL 



The first class of systems to which we shall now apply our general results will be that 

described by the relativistic linear sigma model at finite density (or chemical potential). 

I I 

This was already studied in Ref. [10] where we found that the chemical potential induces 
mixing terms in the Lagrangian and the excitation spectrum can be determined only upon 
an appropriate diagonalization of the matrix propagator. In this section we shall see that 
this complication may be circumvented by the use of Eq. ([3]), at least in the case of type-II 
GBs. We start with a detailed analysis of a simple example in order to illustrate the validity 
of the formulas (TjJ and (Q. 



A. Model with SU(2) x U(l) symmetry 



We recall the mode 
dense quark matter 



that was used to describe kaon condensation in the CFL phase of 
|. It is defined by the Lagrangian, 



C = D^Dy - M 2 0V - A(0Vr, (7) 

where the scalar field <fi transforms as a complex doublet of the global SU(2) symmetry. In 
order to account for finite density, the chemical potential fi associated with the global U(l) 
invariance (particle number) is introduced via the covariant derivative, Do<p = (<% — i/j,)<p. 

Once the value of the chemical potential exceeds the mass parameter M, the scalar 
field condenses. The ground state may be chosen so that only the lower component of <j) 



develops expectation value. At tree level, (02) = v/y/2, where v = a/ — M 2 )/X. This 
breaks the original SU(2) x U(l) symmetry to its U(1)q subgroup, generated by the matrix 
Q = |(1 + T3). The spectrum contains two (as opposed to three broken generators) GBs. 
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The type-II GB is annihilated by 0i and its exact (tree-level) dispersion relation is given by 

Vk 2 + /i 2 -/i. (8) 



10 



The type-I GB is annihilated by a (energy-dependent) linear combination of 02 and 02, and 
its dispersion relation corresponds to the "— " case of 



ojI (k) = k 2 + - M 2 ± Wk 2 + (3/i 2 - M 2 ) 2 . (9) 



After shifting the field 02 by its vacuum expectation value, the Noether currents j% 
associated with the SU(2) symmetry of the model pick up terms linear in 0. From these, 
we may readily express the current transition amplitudes related to the type-II GB state 
\G(k)) as 

(0|jf |G(k)> = -im\G(k)) = + 2^°)(O|0 1 (O)|G(k)). (10) 

Note that this, according to Eq. (j3j), immediately leads to the dispersion relation u = k 2 /2/i, 
which is indeed the low- momentum limit of the exact result jSJ). Moreover, evaluating the 
amplitude (O|0i(O)|G(k)) explicitly gives 

dG* = -2z/iV|(O|0 1 (O)|G(k))| 2 = -ifiv 2 

at |k| = 0. This is in accord with the density rule ([6]) for the ground-state isospin density 
is simply given by the constant term in j®, (0\j^\0) = —fiv 2 , and the relevant structure 
constant of the SU(2) group in the basis of Pauli matrices is /123 = 2. 

Analogously, the coupling of the type-I GB state |7r(k)) to the current is found to be 

<0|j3 M Wk)> = -| [(*" - 2^°)(O|0 2 (O)|7r(k)) - + 2^°)(O|0 2 (O)|7r(k))" . 

The basic complication as opposed to the type-II GB is that now the current amplitude is not 
proportional to a single scalar field amplitude as in Eq. (jlOp . but contains two entangled 
amplitudes. Therefore, these do not simply drop in the ratio as before and have to be 
evaluated explicitly. Using the expressions for the scalar field amplitudes derived in Ref. 



the final result for the transition amplitudes at |k| — > becomes 



(0|j 3 U |7r(k)) = --=V / |k 



2 A J2 \ -1/4 



^2 v ' ' V3/W 2 -M 2 

(0|j s |7T(k)) = - 



/i 2 — M 2 x 1/4 



V2 J\k\ V3/i 2 -M 2 
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A simple ratio of the coefficients of the spatial and temporal current transition amplitudes 
gives the dispersion relation 

/ At 2 - M 2 

which is the correct low- momentum limit of the full dispersion (Q. 
B. General case 

We checked the formulas ([3]) and ([6]) on the simple example of the linear sigma model 
with an SU(2) x U(l) symmetry. We saw that their application is particularly simple for 
the type-II GBs, which carry unbroken charge. On the other hand, the application to type-I 
GBs is obscured by mixing of scalar field transition amplitudes. 

With this experience we are ready to investigate the general linear sigma model with 



arbitrary symmetry. In Ref. 10j, we dealt with the Lagrangian 

C = D^Dy-V(<j>), 

where V{4>) is the most general static potential containing terms up to fourth order in that 
is invariant under the prescribed symmetry. The covariant derivative, -D M = (<9 M — iA^)(j), 
includes the chemical potential(s) assigned to one or more mutually commuting generators 
of the symmetry group. 

Once the chemical potential is large enough, the scalar field develops nonzero expectation 
value, 0o, and must be reparameterized as 

0( x )= ^[00 + ( H ) 

The field U(x) is a linear combination of the broken generators and includes the Goldstone 
degrees of freedom. The "radial" degrees of freedom are described by H(x). 
The Noether current associated with a particular conserved charge T a reads 

f a = -2 Im0 t T a ^0 = -2 Im0 t T a ^0 + 20 f T a ^0. 

Now we just have to insert the parameterization ffTTT) and retain the terms linear in the field 

n, 

Ja,n,Un = -<Pl{Ta, <^II}0O + 2*0 o ^[T a , n]0 O . 



s 



According to Eq. (jlj), this is all we need to calculate the dispersions of the type-II GBs 
present in the system. The commutator and anticommutator will give factors that are purely 
geometrical, i.e., determined by the structure of the symmetry group and the symmetry 
breaking pattern. Without having to evaluate them explicitly, we simply observe that the 
constant (G a ) term is proportional to the chemical potential. This means that the dispersion 
relation always comes out as u = ck 2 //i, where c is some, yet unknown /i-independent 
constant. 

This constant can be determined in a different manner. Recall that (at least at tree level), 
once there is a (at least discrete) symmetry that prevents the potential V(<f>) from picking 
a cubic term, the phase transitions between the normal and the Bose-Einstein condensed 
phases as well as those between two ordered phases are of second order. This means that 
the dispersion relations of the excitations must be continuous across the phase transitions. 
Hence, due to their particularly simple form, they can be matched to the dispersion in the 
normal phase. In the normal phase, a particle carrying charge Q of the symmetry equipped 
with the chemical potential [i, has dispersion to = Vk 2 + M 2 — [lQ, M being its gap at 
[i = 0. At M = [iQ the particle becomes a GB of the spontaneously broken symmetry and 
its low-momentum dispersion relation reads 

k 2 , , 

u= — -. 12 
2[iQ y 1 

Due to the argument sketched above, this dispersion relation then persists to the broken- 
symmetry phase for arbitrary values of the chemical potential. 

IV. NAMBU JONA-LASINIO MODEL 

Now we turn our attention to models of the Nambu-Jona-Lasinio (NJL) type, i.e., purely 
fermionic models with local, nonderivative four-fermion interaction. There, symmetry is 
broken dynamically by a vacuum expectation value of a composite operator. At the same 
time, the GB of the spontaneously broken symmetry is a bound state of the elementary 
fermions, and it manifests itself as a pole in the correlation functions of the appropriate 
composite operators. 

The standard way to calculate the dispersion relation of the GB is as follows (see e.g. 



Ref. [4 1). One first performs the Hubbard-Stratonovich transformation, including the cor- 
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responding pairing channels. The purely scalar action is next expanded to second order in 
the auxiliary scalar field so as to get the GB propagator and find its pole or, equivalently, 
the zero of its inverse. However, the inverse propagator contains a constant term which 
prevents us from identifying the gapless pole directly. This term is removed with the help of 
the gap equation. Only after then may the inverse propagator be expanded to second power 
in momentum in order to establish the dispersion relation. 

The approach proposed in this paper simplifies the calculation of the GB dispersion re- 
lation in three aspects. First, there is no need to use the gap equation in order to establish 
the existence of a gapless excitation. Indeed, since the only assumption, besides rotational 
invariance, is the current conservation, the broken symmetry is already built into the for- 
malism. Second, the loop integrals are only expanded to first order in the momentum, again 
since one power of momentum is already brought about by the current conservation. This 
saves us one order in the Taylor expansion and thus reduces the analytic formulas consider- 
ably. (Of course, the final result must be the same whatever method we use. We just choose 
the less elaborate one.) Third, at least when calculating the dispersion relation of a type-II 
GB, we may take the zero-th component of the external momentum equal to zero from the 
very beginning. This means that, at finite temperature, we get a direct expression of the 
GB dispersion in terms of thermal loops, without having to analytically continue the result 
back to Minkowski spacetime afterwards. 

A. Dispersion relation from fermion loops 

In the models of the NJL type, the Noether currents are fermion bilinears. Their couplings 
to the Goldstone states are therefore given by the one-loop diagrams, symbolically, 



The solid circles denote the full fermion propagators, which are in turn obtained by solving 
the gap equation, while the empty circle stands for the yet unknown vertex coupling the two 
fermions to the GB state. In case of a type-II GB the dispersion relation is then, according 
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to Eqs. and (j4]), given by the ratio, 



W (k) 



k ■ j 



J, 



i'0 



a 




(13) 



where the symbols standing in front of the loops denote the operators inserted in place of 
the crosses. 

Before we explain in detail how the Goldstone-fermion vertex may be extracted from 
another Ward identity, let us note that the explicit evaluation of this vertex may in fact 
often be circumvented. This will be demonstrated later on explicit examples. In general, 
the auxiliary scalar fields introduced by the Hubbard-Stratonovich transformation mix so 
that one has to look for poles of a matrix propagator. However, it seems to be a generic 
feature of type-II GBs that they carry charge of some unbroken symmetry. As a result, they 
may couple to just one of the auxiliary scalar fields and no mixing appears. (For instance, 
recall the model of Sec. IIII Al where the type-II GB couples just to 4>±, while the type-I GB 
couples to both <p2 and <p' 2 -) I n other words, we may devise an interpolating field ip^Xi/j, 
X being some matrix in the flavor space, and calculate instead of the transition amplitude 
(0|j M |7r(k)) the correlation function of and ifi'Xif). The point is that, by the Kallen- 
Lehmann representation, this correlation function is proportional to the desired transition 
amplitude, and that the GB propagator as well as the unknown amplitude (7r(k)\^Xif)\0) 
drop out in the ratio (fl3"j) . 

In practice this means that, instead of the unknown Goldstone-fermion vertex in the loops 
in Eq. (fT3l) . we may simply insert the operator ip^Xip. In particular, in the nonrelativistic 
NJL model, the current j% associated to the charge T a reads 



where m is the common mass of the degenerate fermion flavors. Inserting these operators 
into the fermion loops in Eq. ([TBI , and making the Taylor expansion of the spatial loop to 
the first order in the external momentum, we arrive at the result 



f a =^T a lP, j 



a 



— (^T a Vtlj-V^T a tlj), 




(14) 
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where the trace is performed over the flavor space and S denotes the full fermion propagator. 
At finite temperature l/ft the integral over q° in Eq. (fT4l) becomes a sum over fermion 
Matsubara frequencies, u n = (2n + 1)tt/(3. The GB thus naturally remains gapless up to 
a possible thermal symmetry restoring phase transition, but otherwise its dispersion can 
depend on temperature. 

B. GB fermion coupling from Ward identity 

As already stressed in the discussion below Eq. ffT3l) . the formula (fT4|) has a limited use. 
It applies only to the cases in which the GB couples to just one of the auxiliary scalar fields 
introduced by the Hubbard-Stratonovich transformation (no mixing). Unless this happens, 
the (unknown) couplings of the GB to the fermion bilinear operators do not drop in the 
ratio (TT31). 

We also noted that this no-mixing property is most likely to be generic for type-II GBs, 
which carry charge of some unbroken symmetry. Nevertheless, in any case, the GB-fermion 
coupling can be evaluated directly by using yet another Ward identity [ijj]. Let us introduce 
the connected vertex function, G£(x,y,z) = (0\T{j£(x)ip(y)ip^(z)}\0). Taking the diver- 
gence and stripping off the full fermion propagators corresponding to the external legs, we 
find the Ward identity for the proper vertex, 



This vertex exhibits a pole due to the propagation of the intermediate GB state. Near 
the GB mass shell, it may be approximated by a sum of the bare vertex and the pole 
contribution [ijj]. Writing the inverse nonrelativistic fermion propagator as 5 ,_1 (u;,p) = 

2 

lu — §^ + ^ — S(u;,p), where E is the self-energy, the bare part of the propagator exactly 
cancels against the bare vertex in the Ward identity ( |T5l) . The remaining term expresses 
the residuum at the GB pole in the vertex function — which is proportional to the desired 
GB-fermion coupling — in terms of the fermion self-energy, so that 



The normalization constant N is proportional to the coupling of the GB to the broken 
current itself, but this is not important since this factor only appears in the ratio ffT3l) . 



k^{p + k,p) = T a S~\p + k) - S-\p)T a . 



(15) 




(16) 
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The GB-fermion coupling is apparently determined by the symmetry-breaking part of the 
fermion self-energy. Once this is calculated, e.g., in the mean-field approximation, it may 
be plugged into the loops in Eq. ([TBI) to determine the GB dispersion. 

In conclusion of our general discussion let us remark that all the results presented so far 
in Sec. IIVI extend straightforwardly to cases including Cooper pairing of fermions. Once 
the GB is annihilated by a difermion operator instead of a fermion-antifermion one, we just 
have to flip the direction of some of the arrows in the fermion loops. The S(q) in Eq. (TH]) 
then denotes the matrix propagator in both the flavor and the Nambu (ip, space. 

In addition, Eqs. ( lT5i) and ( fl6|) modify accordingly: The matrix T a must be replaced 
with the block matrix in the Nambu space, ( a _ t t j ■ In the end, the GB-fermion coupling 
vertex comes out as 

=^ = [^4p)tI + r.s^(p + k)] , (17) 

where denotes the anomalous (pairing) part of the fermion self-energy. 



C. Three-flavor degenerate Fermi gas 

The nonrelativistic Fermi gas consisting of three degenerate fermion species has been 
investigated theoretically in great detail 



16l . Il7| . The reason is that such a system has 
already been successfully realized in experiments using atoms of 6 Li and 40 K. At the same 
time it provides an interesting analogy with quantum chromodynamics (QCD). A recent 
theoretical analysis suggests that under certain conditions, bound clusters of three atoms 
may form, being; singlets with respect to the global SU(3) symmetry, very much like baryons 
in QCD [l8(. Such a three-flavor atomic Fermi gas thus creates a natural playground for 
testing properties of quark matter. 

In the following, we are going to illustrate the utility of our general formulas and rederive 
the results of Ref. js| concerning the type-II GBs in the three-flavor system. To that end, 
recall that with a proper choice of the interaction, the global SU(3) symmetry is broken by 
a difermion condensate, (ipify). Since this is antisymmetric under the exchange of it 
transforms in the 3 representation of SU(3). This leads to the symmetry breaking pattern 
SU(3) x U(l) -> SU(2) x U(l). 

Once we choose the condensate to point in the third direction, the generator A 8 acquires 

13 



FIG. 1: Fermion loop for the calculation of the type-II GB dispersion in the three- flavor Fermi 
gas. The numbers adjacent to the fermion propagators denote the respective flavors. The upper 
line represents the normal propagator of the (ungapped) third fermion flavor, while the lower line 
represents the pairing of the first two flavors. 

nonzero vacuum expectation value. Based on the general discussion in Sec. [TTJ, we anticipate 
that the four broken generators, A4, A5, A6, A7, couple to a doublet of type-II GBs with 
quadratic dispersion relations. Due to the unbroken symmetry, these two GBs have the 
same dispersion. For the explicit calculation, we choose the GB in the sector (A 4 , A 5 ). 

To determine its dispersion relation, we need to evaluate the loop in Fig. [1] at nonzero 
external momentum k. [In order to make contact with literature [8(, we do not perform 
the Taylor expansion in the external momentum as in Eq. (1T4"1) .] Note that since the type- 
II GB couples to two broken currents, j% and j'5 , we are free to choose any of them to 
calculate its dispersion. Here, we take up j% . The two propagators that appear in the loop 
are given by the finite-temperature expressions (see Ref. [8j) ^ 33 (2z/ n ,k) = \j{iv n — £ k ), 
J r 12 (i^ n , k) = — zA/[(iz/ n ) 2 — E%\, where A is the gap parameter, £k = k 2 /2m — //is the bare 
fermion energy measured from the Fermi surface, and = y^jf-KA 2 . The interpolating 
field for the GB is chosen as ip T \-jip. This may be understood as follows. By our choice, the 
ground state condensate carries the third (anti)color of the representation 3. The current 
j% whose transition amplitude we calculate, then excites a difermion state with the first 
(anti) color. The GB thus carries the quantum numbers of a pair of fermions with the 
second and the third color, respectively. The same conclusion follows immediately from Eq. 
(JTTj): The fermion self-energy is proportional to A2 and {A2, A4} = — A7. 

The point of using Eq. 0131) to calculate the GB dispersion relation is that all common 
algebraic factors including the trace over the flavor space cancel in the ratio. After stripping 
off such trivial factors, we are thus left with the loop integral 




(^n) 2 - m 



1 



q-k 




1 



(18) 



14 



where the two options in the curly brackets stand for the insertion of the temporal and spatial 
components of the current, respectively. Note that the Matsubara frequency entering the 
loop graphs has been set to zero because, as already remarked above, in case of a type-II 
GB the transition amplitude of j° is nonzero even at vanishing energy. 
After performing the Matsubara sum, Eq. (1181) acquires the form 



J(k) 



d 3 q 



-k+2q 

2 m 



IE, 



q-k 



f(E t 



q-kj 



E, 



q-k 



HJ + /(-^q-k) 



E, 



q-k 



£q 



(19) 



where f(x) = l/(e^ x +l) is the Fermi-Dirac distribution. This result is completely equivalent 
to Eq. (A9) in Ref. jg], though not identical, of course, for we calculate a different loop 
graph here. The GB dispersion may be directly inferred from here by means of Eq. (fl3l) . 
We omit the explicit expression since it is rather cumbersome and not particularly revealing. 



V. FOUR-FLAVOR DEGENERATE FERMI GAS 



In this section we are going to discuss the nonrelativistic degenerate Fermi gas with 
four fermion flavors which, to our best knowledge, has not been analyzed in literature so 
far. Besides being the next-to-simplest nontrivial case of a nonrelativistic fermionic system 
exhibiting type-II GBs, it als o provides an interesting analogy with the two-color QCD with 



two quark flavors 
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2i 



We start by assuming the existence of four degenerate fermion species with a common 
chemical potential /i, which interact by a contact four-fermion interaction invariant under 
the global SU(4) x U(l) symmetry. For the time being, all we need to know is that this 
interaction destabilizes the Fermi sea by the formation of Cooper pairs. The detailed form 
of the interaction will be specified later. 



A. Symmetry breaking patterns 

First we analyze the symmetry properties of the ordered ground state. Analogously to 
the three-flavor case, the original SU(4) x U(l) symmetry is spontaneously broken by a 
difermion condensate (ipiipj), which transforms as an antisymmetric rank-2 tensor of SU(4). 

The analysis of the symmetry breaking patterns is much simplified by noting the Lie 
algebra isomorphism SU(4) ~ SO (6). The six- dimensional antisymmetric tensor represen- 
tation of SU(4) is equivalent to the vector representation of SO(6). This correspondence is 



15 



worked out in detail in the Appendix. The order parameter is thus represented by a complex 
6- vector, 0, or equivalently, by two real 6- vectors, <pR = Re0 and 0/ = Im0. The situation 
is pretty much similar to that in spin-one color superconductors where the order parameter 



is a complex 3- vector of SO (3) 22j, |23|. \24\. 



There are two qualitatively different possibilities, distinguished by the relative orien- 
tation of <pR and 0/. First, if 0r and 0/ are parallel, there is obviously an unbroken 
SO(5) rotational invariance. The symmetry breaking pattern is SO(6) x U(l) — > SO(5), 
or SU(4) x U(l) — ■> Sp(4), with six broken generators, transforming as a 5-plet and a sin- 
glet of the unbroken SO(5), respectively. Note that in this case, the U(l) symmetry of 
the Lagrangian may be used to even eliminate entirely the imaginary part of 0, i.e., to set 
<ft i = 0. Let us remark that it is this phase of the nonrelativistic four-flavor Fermi gas that is 
analogous to the two- flavor two-color QCD. There, the reality of is enforced by C and P 
conservation. Also, there is no global U(l) symmetry in two-color QCD as the corresponding 
phase transformations are the axial ones, being broken by the axial anomaly 20, l2lj ]. 

The more interesting, but also more complicated possibility is that with nonzero 0# and 
0/ pointing in different directions. Then, the SO(6) x U(l) symmetry may be exploited 
to cast the order parameter in the form = (a,j3 + ry, 0, 0, 0, 0) T , with real a, (3, 7, and 
both a and 7 nonzero. In general the SO(6) x U(l) symmetry is now broken down to 
SO (4), thus breaking 10 generators, forming two vectors of the unbroken SO (4) as well as 
two singlets. However, when (3 = and |a| = I7I, there is an additional unbroken U(l) 
symmetry generated by a combination of a rotation in the (<Pr, 4>i) plane and the original 
U(l) phase transformation. There is thus one less broken generator and one less GB. 



B. Fermion propagator and excitation dispersions 

According to our general discussion above, we need to know the charge densities in the 
ground state in order to determine the spectrum of GBs of the broken symmetry. The charge 
densities, (ip^Taip), are obtained from the trace of the fermion propagator with the charge 
matrix T a . At this place, we have to specify the Lagrangian of the system. 

Anticipating the Cooper pairing of the fermions, we choose the interaction Lagrangian in 
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a form suitable for a mean-field calculation, 

C = ^ (id t + + fi) V + G^2 iV^Vf, (20) 

where if) is the four-component fermion field and Ej are a set of six basic antisymmetric 4x4 
matrices, transforming as a vector of SO (6), see Eqs. 11 All) and (]A2j) . Thus, the Lagrangian 
(120]) is manifestly SO(6) ~ SU(4) invariant. 

Upon introduction of an auxiliary scalar field fa by the Hubbard-Stratonovich transfor- 
mation, the Lagrangian becomes 

C = tf(id t + ^ + f ?)l>- -^4>]fa + + H.c. (21) 

The vacuum expectation value of provides the order parameter for symmetry breaking. 
In the mean-field approximation, is simply replaced with its vacuum expectation value 
and the semi-bosonized Lagrangian 112 lj) then yields the inverse fermion propagator in the 
Nambu space, 

$ + (Go)' 

where Gq are the bare propagators and $~ = faTi\, $ + = ( < 1 )_ )^. Inverting this block 2x2 



G- 1 



matrix, we obtain 




G = 
where 

G ± = [(G^y 1 - $ T G F <1> ± ] _1 , E ± = -G^ ± G ± . (22) 

At finite temperature, the bare thermal propagators are simply Gq {iv n ^ k) = 1/ (iv n — £ k ) 
and Gq (iu n ,k) = \j{iv n + £ k ). As explained in the previous section, we can exploit the 
SO(6) symmetry of the Lagrangian fl20]) so that only two components of the 6-vector <p are 
actually nonzero. Let us denote them generally as a,b (a ^ b), so that $~ = a S^ + faY*\. 
In the previous section we simply took a = 1 and 6 = 2, but when later calculating the 
charge densities, it will be more convenient to set a = 1 and 6 = 6. 

According to Eq. (|22|) . the full mean-field fermion propagator G~ reads 

G~(iu n , k) = -K - £ k )(^ + £ + $+$-)- 1 . (23) 
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Using the properties of the matrices we have 

= \<p a \ 2 + \<p b \ 2 - 2iIm(0 a <^)E a Ej = \<p a \ 2 + \<p b \ 2 - 2Im(0 a 0*)[P a ( 6 +) - P^]. 

The projectors pf^ are defined in Eq. (1A5[) . With projectors in the denominator, the 
matrix in the expression (1231) for G~ can easily be inverted and yields the final result 

G~(iv n ,k) = -{iv n - f k ) 





2 +l< 




2 + 2Im(0 a 0*)[P a ( + ) -P a ( - ) ] 


+ £ + 1*. 


2 +106 


2)2-4[Im(0 a 0*P 



p(+) p(-) 

aft _|_ ab 



^ + + \<Pa\ 2 + N 2 - 2 lm(0 a 0*) + « + \<Pa\ 2 + |<M 2 + 2 lm(0 a 0*) 

(24) 

From Eq. (|24p we can see that there are two doublets of fermionic excitations over the 
ordered ground state, with generally nonzero gaps. Their dispersion relations are most 
conveniently written in the form 



J ± \k) = ^/a + \<PaTtfo\ 2 . (25) 

Two special cases deserve particular attention. First, in the SO (5) symmetric phase is set 
equal to zero so that all four excitations are degenerate, forming a fundamental quadruplet 
of SO(5) ~ Sp(4). Second, in the SO(4) x U(l) symmetric phase, \<f> a \ = |</>&| with their 
complex phases differing by tt/2. Eq. (T25]) then implies that two of the fermions are actually 
gapless. Altogether we have two pairs of degenerate fermions, which transform as doublets 
under the respective SU(2) factors of SO(4) ~ SU(2) x SU(2). 



C. Charge densities and GBs 

The ground state densities of the conserved charges are obtained from the propagator 

G~ as 

d 3 k 



For sake of this section, it is convenient to fix the direction of the order parameter <j) in the 
SO(6) space. As already remarked above, we choose a = 1 and 6 = 6. With the explicit 
form of the Ej matrices given in Eq. (1A2I) . we find 

?Ti yt _ P (+) _ p(-) _ _ 
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Since ^SiSg = Mi 6 is one of the SO (6) generators My defined in Eq. (|A6|) . the orthogonality 
relation (1A8I) implies that it can also be the only one which acquires nonzero density in the 
ground state, 

(M 16 ) = 8 £ / ^ + £2 + |^|2% 6 |^ _ 4[imT«f 



d 3 k , 



l-2/M+)(k)) l-2/(o;H(k)) 
w(+)(k) wH(k) 



(2tt) 3 ~ 

after the Matsubara summation. 

The spectrum of GBs may be understood as follows. First, in the SO(5) symmetric phase, 
6 = and there is no charge density in the ground state. According to our general results, 
there are then six type-I GBs — five corresponding to the coset SO(6)/SO(5), and one to the 
broken U(l). 

From now on, we shall analyze the SO (4) symmetric phase of the model ( 120]) . i.e., assume 
that both 0! and 06 are nonzero and have different complex phases. The condensate 0i 
breaks all SO(6) generators of the form Mu, where i = 2, ... ,6. On the other hand, the 
condensate 06 breaks all SO(6) generators of the form Mjq, where j = 1,...,5. This is 
altogether nine spontaneously broken generators of SO (6). One of them, M 16 , is a singlet 
of the unbroken SO (4). The remaining eight ones group into two SO (4) vectors, Mu and 
Mj6, with i — 2, ... ,5. Note that, by Eq. (1A7I) . [Mu, Mj 6 ] = 2iMi 6 (no summation over i!). 
Since M 16 has nonzero density in the ground state, the two broken generators Mu and M ie 
couple to a type-II GB with a quadratic dispersion relation. These eight broken generators 
are thus associated to an SO (4) vector of type-II GBs. 

Finally, for 6 = ±i0i there is a residual unbroken U(l) symmetry generated by 1 ±M 16 . 
Note that the combinations Mu ± iM^ are charged under this unbroken symmetry, carrying 
opposite charges. It may also be illuminating to observe that, for instance, when 06 = —i<f>i, 
the matrix $~ has the form 

/ 1 o\ 

-10 
\ / 

In words, only the first and third fermion flavors participate in the pairing. The second and 
fourth do not, and represent the two gapless excitations predicted by Eq. (1251) . 



$~ = 20! 
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G-(q) 




FIG. 2: Fermion loop for the calculation of the type-II GB dispersion in the four-flavor Fermi gas 
described by the Lagrangian ([20]) . 



Let us now calculate the dispersion relation of the type-II GBs in this SO(4) x U(l) 
symmetric phase. (As will be shown in the next section, the relation 6 = ±z0i is one of the 
two possible solutions of the mean-field gap equation, the other one corresponding to the 
SO(5) symmetric phase.) As already noted above, there are four type-II GBs, that couple 
to the broken generators Mu and Since these form a multiplet of the unbroken SO(4), 
they have all the same dispersions. We choose i — 3 for which the broken generators have a 
particularly simple form, 



Mi 3 = - , M 36 





Together with M 16 these generate an SU(2) subalgebra of SU(4). 

The generators M i3 and M 36 must both couple to the sought type-II GB. When acting 
on the ground state, they both excite the third component of 0, so the interpolating field for 
the GB will be 4 ,T ^4 > - 111 order to determine the GB dispersion, we evaluate its coupling 
to the broken current of M 13 , which is given by the loop graph in Fig. [2l analogous to that 
for the thee-flavor case in Fig. [TJ After performing the Matsubara sum, and throwing away 
trivial numerical factors, we arrive at the loop integral, analogous to Eq. (|19|) . 



J(k) 



d 3 q r ii e 



,,(±) , ,,(T) ,,(±) ( .(T) 



where the upper and lower signs were obtained with 6 = respectively. Ultimately, 

both yield the same dispersion, which is found by the Taylor expansion to first order in 
momentum and plugging the result into Eq. (jlj). 
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D. Gap equation 



The actual value of the order parameter is determined by the gap equation. To that 
end, we need to know the anomalous part of the matrix propagator which follows from Eq. 

E~(iu n , k) = -- — — — (0 a Sj + fe £j)G - O n , k). 

In the mean-field approximation, we require the self-consistency at the one-loop level, which 
leads to the gap equation 

d 3 k 



1 f MAT 

^ = 2G pJ2J ^tTr[S 4 -(^,k)]. 



After working out the algebraic trace and the Matsubara sum, this set of gap equations for 
cf) a and 4>b may be rewritten in a particularly symmetric form, 

<f> a + i<f> b = 4G{<f> a + i<l> b )I^ {</>), 

(26) 

a -^ = 4G(0 a -# 6 )/ (+) (0), 

where 

W) J (2tt) 3 w (±)(k) ' 
The gap equations (126]) apparently admit two qualitatively different solutions. The first 

one corresponds to lm(0 a 0£) = 0, i.e., the SO(5) symmetric phase. The gap equation then 
simply reads 

1 =4G/ (±) (0). (27) 

Either sign can be used in the integrals since, in this case, the o/^ dispersions are degenerate. 

The second solution has lm(0 a 0£) ^ 0. Eq. ( 1261) then implies that <p a = ±«06- For the 
two possibilities, the gap equations are 

l=4G7 ( - ) (0), for0 a = +# b , 

(28) 

1 = 4GJ«(0), for0 a = -i0 6 . 

Note that this solution of the gap equations (I2"6"j) precisely corresponds to the 
SO(4) x U(l) phase discussed in Sec. IV Al In either case, <fi a = ±z0 b , only the two gapped 
fermion flavors contribute to the gap equation (1281) . 

It is interesting to observe that the mean-field analysis does not support the most general 
form of the order parameter <fi. Rather, either the real and imaginary parts of are parallel, 
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breaking the SO (6) x U(l) symmetry to SO (5), or they are perpendicular with the same 
magnitude, breaking the symmetry to S0(4) x U(l). 

Of course, which of these two possibilities is the actual ground state may depend on the 
magnitude of the coupling constant and the chemical potential, and can only be decided 
by solving the gap equations (1271) and (1281) and comparing their free energies. This task 
is however, beyond the scope of the present paper where the model (|20|) only serves as an 
illustration of the general features of symmetry breaking in presence of charge density 

VI. PSEUDO-GOLDSTONE BOSONS 

So far we have been dealing exclusively with exact symmetries for which, once sponta- 
neously broken, the Goldstone theorem ensures the existence of exactly gapless excitations. 
Real world is, however, more complicated. Almost as a rule, global internal symmetries are 
usually only approximate, leading to approximately conserved currents and, when sponta- 
neously broken, to the existence of excitations with a small gap, so called pseudo-Goldstone 
bosons (pGBs). (In the following, we shall keep using the term Goldstone boson unless the 
existence of a nonzero gap is important for our discussion.) 

Let us just note that, in fact, the only generic type of an exact global symmetry is the 
particle number conservation. This is, however, only Abelian and hence does not lead to 
the nontrivial phenomena studied in this paper. In addition, gauge symmetries are always 
exact. Here, once the gauge is appropriately fixed and the remaining global symmetry 
spontaneously broken, there are no physical GB states as they are "eaten" by the gauge 
bosons, making them massive by the Higgs mechanism. It turns out that the number of 
massive gauge bosons is always equal to the number of broken generators even if the number 
of would-be GBs is smaller [251 ]. 

A. Counting of pseudo-Goldstone bosons 

It should be stressed that for approximate symmetries, there is no rigorous theorem anal- 
ogous to the Goldstone's, which would guarantee the existence of "approximately gapless" 
excitations. On a heuristic level, however, the conclusions of the theorems of Goldstone and 
of Nielsen and Chadha may be carried over to this less perfect, yet more common and physi- 
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cal situation. Before a more detailed discussion of approximate symmetries, let us illustrate 
this on a simple example. 

Both Goldstone and Nielsen-Chadha theorems make predictions about the low- 
momentum limit of the GB dispersion relation without actually specifying what the low 
momentum means. What is the type and number of pGBs associated with a spontaneously 
broken approximate symmetry, therefore depends on the scale at which the symmetry is 
broken explicitly. Consider the simple model studied in Sec. IIII A[ The type-II GB has a 
massive partner which is annihilated by 4>\, i.e., its antiparticle, and their exact tree- level 
dispersions may be jointly expressed as 



see Eq. (jSJ). At low momentum, we clearly see a gapless mode with a quadratic dispersion 
and a mode with the gap 2/x. On the other hand, at high momentum both modes have 
almost linear dispersion. 

So, once the SU(2) x U(l) symmetry of the model (0) is explicitly broken, the resulting 
pGB spectrum depends on whether the scale of the explicit breaking is smaller or larger 
than the chemical potential fi. For weak breaking, we still see an approximately gapless 
mode with a quadratic dispersion and a mode with a gap approximately 2/i. On the other 
hand, for strong explicit breaking, larger than /i, we can no longer resolve the quadratic 
onset of the dispersion, and instead find two pGBs with linear dispersion relations. Note 
that in both cases the Nielsen-Chadha counting rule is satisfied even for pGBs. 

B. Ward identity and dispersion relations 

Let us now consider a symmetry which is explicitly broken by adding a small term to 
the Lagrangian of the theory. The associated current then becomes only approximately 
conserved and we find d^j^ = A a C, where A a C is the variation of the Lagrangian under 
the symmetry transformation. The basic assumption of an "approximate" spontaneously 
broken symmetry is that there is still a mode |vr(k)) which couples to the current j%. The 
parameterization (j2J), which is based just on rotational invariance, thus still holds. What 
does change is, of course, Eq. (jHJ) which now becomes 



u = 



>/k 2 + fi 2 ± fi, 




(29) 
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Up to this change, we proceed along the same line of argument as in Sec. [TTJ The Ward 
identity (j3J) acquires a new term that accounts for the explicit symmetry breaking, 

d;(0\T{f a (x)f b (y)}\0) = (0\T{A a £(x)f b (yW) + ~ V^MWMW- 

Since both the broken currents and A a C couple to the GB state, both correlators can be 
expressed via the Kallen-Lehmann representation. We find that cancelation of the one- 
particle poles exactly requires the energy and momentum to fulfill the dispersion relation 
(|2"5|) . In the limit |k| — > 0, the Ward identity leads to the modified density rule, 

lm[G a G* b + Lo(0)G a F b *} = i/ afcc (0| Jc °|0), 

cf. Eq. (EJ). The explicit breaking of the symmetry enters here in terms of the (nonzero) 
energy gap, u(0). 

Eq. (|2"§I) may be used to give a precise meaning to the otherwise vague term "approx- 
imately gapless mode", as well as to distinguish type-I and type-II GBs in the case of an 
explicitly broken symmetry. Let us thus, again, assume that G a (|k|) has a nonzero limit as 
|k| — > 0, and denote (0| A a £|7r(k)) = if a (|k|), for sake of brevity. 

According to Eq. (l29il . the energy gap u(0) of the pGB is given by the solution of the 
equation 

Lu 2 (0)F a + u(0)G a = H a , 
which has exactly one positive solution provided all coefficients, F a , G a , H a , are positive, 

"(0) = (~G a + y/G* + 4F a H a ) . 

The low-momentum behavior of the GB dispersion is apparently determined by the rel- 
ative value of G\ and F a H a . If the explicit symmetry breaking is weak enough, i.e., if 
F a H a G 2 al the term quadratic in u may be neglected in comparison with that linear in lo, 
and the low-momentum dispersion relation of the pGB reads 

-(k) = ^ + k 2 ^, (30) 

as opposed to Eq. (TJJ. We thus find a type-II pGB. 

If, on the other hand, F a H a ^> G 2 , the uj 2 term dominates over the uj term in Eq. ( 1291) . 
and the resulting low-momentum dispersion relation, 

u(k) = ^ + k2, (31) 
is one of a type-I pGB with nonzero mass squared H a / F a . 
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C. Nambu Jona-Lasinio model 

The formulas (l29]) - (}3Tj) are absolutely general but, unfortunately, not of much use unless 
we are somehow able to calculate the new amplitude, (0|A a £|7r(k)). Nevertheless, there 
is a wide class of systems where this can be determined in the same way as the current 
amplitudes F a and G a : They are the models of the NJL type where the symmetry-breaking 
operator is bilinear in the fermion fields. 

First of all, note that this is all we need, for instance, to provide a realistic description 
of dense quark matter since both quark masses and the various chemical potentials can be 
taken into account by adding appropriate bilinear operators into the Lagrangian. 

The point is that in such a case, the dispersion relation of the GB f[5Uj) [or fpJTj) as well] 
is given by the ratio of fermion loops, just like in Eq. ffT3l) . which now becomes 




All that was said in Sec. II V Al then holds without change also in the presence of explicit 
symmetry breaking. 

VII. CONCLUSIONS 

In this paper, we analyzed the general problem of spontaneous symmetry breaking in 
presence of nonzero charge density. Based on a Ward identity for the broken symmetry, we 
showed in Sec. [Ill that nonzero density of a commutator of two broken currents implies the 
existence of a GB with nonlinear dispersion relation. Under fairly general circumstances, this 
dispersion is quadratic. The only assumption is the existence of a nonzero low-momentum 
limit of the transition amplitude of the spatial current, F a . 

It is interesting to compare this result to that of Leutwyler |9[]. In his approach, nonzero 
density of a non-Abelian charge leads to the presence of a term with a single time derivative 
in the effective Lagrangian. Together with the standard two-derivative kinetic term, this 
results in a quadratic dispersion of the GB. In fact, however, Leutwyler's kinetic term is 
proportional to the same constant F a that we use here just to parameterize the current 

25 



transition amplitude. Thus, regarding the possibility of the existence of GBs with energy 
proportional to some higher power of momentum, our results are equivalent. Nevertheless, 
the comparison with Leutwyler's approach gives us new insight into this open question: In 
order to have a GB with the power of momentum higher than two in the dispersion relation, 
the standard kinetic term in the low-energy effective Lagrangian would have to vanish. It 
is hard to imagine how such a theory would be defined even perturbatively, yet such a 
possibility is not theoretically ruled out. 

In Sec. II I II we demonstrated the general results of Sec. [Ill on the linear sigma model, 
which we studied in detail in our previous work [lo|, . We were thus able to prove on a 
very general ground, that the tree-level dispersion of type-II GBs in such a class of models 
is of the form (fT2"j) . 

While in the framework of the linear sigma model we merely verified results that could 
be obtained in a different manner, in the models of the NJL type our general results are 
practically useful. In Sec. [TV] we suggested to calculate the GB dispersion relation by a 
direct evaluation of its coupling to the broken current. This procedure considerably reduces 
both the length and the complexity of the analytic computation as compared with the 
conventional approach using the Hubbard-Stratonovich transformation and the consecutive 
diagonalization of the quadratic part of the action. 

In Sec. |V]we analyzed a particular example of a system exhibiting type-II GBs: The non- 
relativistic gas of four fermion flavors with an SU(4) x U(l) global symmetry. We showed 
that the mean- field gap equation admits two qualitatively different solutions, i.e., two dif- 
ferent ordered phases. In the first one, all fermion flavors pair in a symmetric way, leaving 
unbroken an SO (5) subgroup of the original symmetry. The six broken generators give rise 
to six type-I GBs. 

In the other phase, only two of the four fermion flavors participate in the pairing. The 
other two flavors remain gapless. In this case, the unbroken symmetry is SO(4) x U(l) 
and one of the broken charges — corresponding to the difference of the number of the paired 
and unpaired fermions — acquires nonzero density. Consequently, eight of the nine broken 
generators give rise to four type-II GBs with a quadratic dispersion. Which of these two 
phases is the actual ground state of the system, can only be decided upon the solution of 
the gap equation. 

Finally, in Sec. I VII we extended the results achieved in the preceding parts of the paper to 
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spontaneously broken approximate symmetries. We discovered that once the notions of the 
type and low-momentum dispersion relation of the GB are defined and treated carefully, the 
general results remain valid. In particular, in the NJL model this allows their application 
to the realistic description of dense quark matter. 
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APPENDIX A: ALGEBRA OF ANTISYMMETRIC 4x4 MATRICES 



In Ref. [19] we showed that a general unitary antisymmetric 4x4 matrix E with unit 
determinant may be cast in the form £ = rijEj, % — 1, . . . , 6, where £j is a set of six basis 
matrices, satisfying the constraint 

eJe,- + £t£, = 26ijl, (Al) 

and rti is an either real or pure imaginary unit vector. By this construction we provided 
an explicit mapping between the cosets SU(4)/Sp(4) and SO(6)/SO(5). We also found a 
particular explicit solution to Eq. (1A1I) . 

r 2 

(A2) 

%T 2 
%T 2 

By a simple complex conjugation, using the antisymmetry of the E^s, Eq. (I A II) translates 
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to 



E^Ej + Ej-E| = 25«1. (A3) 
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In the context of the present paper, n« no is longer a unit vector, but rather an arbitrary 
complex six- dimensional vector. We start by showing that its real and imaginary parts 
realize independent real representations of the group SU(4) ~ SO (6). This acts on E as 
E — > UHU T . For an infinitesimal transformation, U — 1 + iX with a traceless, Hermitian 
X. Using the orthogonality relation 

Tr EjEj = 45ij, (A4) 

we derive the change of the coordinate rij induced by the transformation U, 

8m = ^Tr[Sl(XE, + E,X^)] = ^Tr^Ef). 

(We transposed the second term and used the antisymmetry of Ej.) Now for traceless X, 
the matrix iTr(XEjE|) is, with a little help from Eq. (1A3I) . readily seen to be real. As an 
immediate consequence, the real part of rij transforms into the real part and analogously for 
the imaginary part, as was to be proven. 

In fact, just the six matrices Ej generate the whole algebra of SU(4). To see this, note 
that iSjEj, i 7^ j, are 15 independent matrices with the following properties: 1. They are 
Hermitian. 2. Their trace is zero. 3. Their square is equal to unit matrix. 

Proof. Hermiticity follows from (iEj-Ej)^ = — iSjEj = iEjEj, by Eq. (IA3I) . The same 
identity together with the unitarity of Ej implies (zE,Ej) 2 = — EjEjEjE^ = EjEjEjEj = 1. 
The tracelessness is an immediate consequence of Eq. ( 1A4I) . 

The fact that the matrices iEjEj square to one tells us that their eigenvalues are ±1 and, 
moreover, both with the twofold degeneracy, in order to ensure zero trace. We may then 
write iEjEj = pjj — \ where the projectors P^- are conveniently expressed as 

4 +) = (S,-^-)(S,-^,) t , 

(A5) 

Pp = (E i + iE i )(E i + zE i ) t . 

The identity ( 1A3I) implies that the product iEjEj is, for i ^ j, antisymmetric under the 
exchange of % and j. To make this antisymmetry explicit, let us define a new set of matrices, 

M iS = |(E,Et - E.-Ej) = i(E,Et - ^1). (A6) 

By means of Eqs. (lAlj) and (1A3|) . these may be shown to satisfy the commutation relation 

[Mij , M kl ) = 2i(5 u M jk + 6 jk M u - 5 ik M 3l - 8 fl M ik ). (A7) 
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The matrices My thus provide an explicit mapping between the Lie algebras of SU(4) and 
S0(6). 

Finally, Eqs. fjAlH . (IA30 . and 0A40 imply the identity, analogous to that for the Dirac 7 
matrices, 

Tr(SjE]E fc E|) = 4(^5^ - S ik 8ji + S a S jk ). 
As a consequence we have the orthogonality relation for the SO (6) generators M^-, 

TrCMyMfeO = 4(<fofy - (A8) 

Note also that none of the identities proven in this Appendix relies on the explicit rep- 
resentation of the matrices Eq. (1A2j) . All the results thus hold independently of the 
particular realization of £j. In a sense, the construction worked out here is similar to that 



of the Clifford algebra of the spinor representations of t 
case of SO(2n), this is realized by matrices of rank 2 n 
matrices for the group SO (6). 



re orthogonal groups. However, in 



261 ] . while here we need just 4x4 
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